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SMEARING FORMULAS FOR DENSITY MATRICES 
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Arnimallee 14, 14195 Berlin, Germany 

We report on the development of a systematic variational perturbation theory for 
the euclidean path integral representation of the density matrix based on new 
smearing formulas for harmonic correlation functions. As a first application, we 
present the lowest-order approximation for the radial distribution function of an 
electron in a hydrogen atom. 
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QJ ' 1 Introduction 



Quantum statistical free energies can be calcul ated to any desired accuracy 
with the help of variational perturbation theory Uddi Th is is a systematic gen- 
"^ ' eralization of an earlier rough variational approach ^^ based on the Jensen- 

Peierls inequality. The approximation consists of an optimized expansion, in 
which each term can be expressed via a simple smearing formula applied to 
powers of the interaction, which can be polynomial as well as nonpolynomial. 
In Sect. [2] we sketch the extension of this theory to density matrices, 
thus enabling us to calculate very accurately local statistical properties of a 
Q ' quantum mechanical system at all temperatures and coupling strengths. The 

O . corresponding smearing formulas are derived in Sect.[3]and applied in Sect. [4] 

to obtain the lowest-order approximation for the density distribution of a 
Coulomb system. 
> 

On 

2 Variational Perturbation Theory for Density Matrices 
(N 

^— ; , Variational perturbation expansions for the free energy approximate arbi- 

trary quantum statistical systems by optimized local perturbation expansions 
around harmonic systems with different trial frequencies. The optimization is 
performed separately for each path average xo = J x(t)/H/3. This ensures a 
rapid convergence of such expansions at higher temperatures by removing the 
fluctuations of xq which diverge linearly with the temperature. The fluctua- 
tions of Xq are accounted for at the end by an ordinary integral. Since these 
^ ' fluctuations probe, at higher temperatures, the entire potential, they must be 

done numerically, and the fact that this is always possible with high accuracy 
is one of the reasons for the quality of the approximations ^. 
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When applying variational perturbation theory to density matrices, the 
special role of the xo-fluctuations diappears since xq remains always close to 
some average of the fixed endpoints of the paths. Then xq may be treated 
perturbatively ^ together with the other Fourier components of the paths. 
We may therefore directly look for an optimized perturbation expansion of 
the path integral for the unnormalized density matrix 

Xb,h/3 

p(x a ,x b ) = / T>x exp< —rA\x] > , (1) 

X a ,0 

where p = 1/kgT and A is the euclidean action. In order to obtain a vari- 
ational approximation, we divide the full action into a harmonic trial action 
with center x m i n — x m ; n (x a , Xb) and frequency fl = Q(x a ,Xb, x m i n ). De- 
noting the trial action by -4n,a; min , the remainder _4i nt = A — Aa tXtnin = 
J dr Vint (x(t)) is treated as a perturbation. The result can be written as 
an exponential of a cumulant expansion (cumulants being indicated by sub- 
script c) 

p~{x a ,X b ) = pa,Xr*ta{%a,Xb) 

x exp j-i ( A„tN }%Z n , c + 7^2 < AM >n"Z.- -■•]• ( 2 ) 

The prefactor /5n, KmJn (a; , Xb) is the unnormalized density matrix of the dis- 
placed euclidean harmonic propagator. The connected correlation functions 
in the exponent of ([2]) consist of the harmonic expectation values 

[A\ u „ .. 

expj --^n, Xmin N f- ■ (3) 

Truncating the series @ after the A^th term, we find the approximation 
PN(x a ,Xb', fl,x nlul ) to the euclidean propagator (TT]). As the exact propaga- 
tor p{x ai Xb) does not depend on the variational parameters, we expect the 
best approximation ppf(x a ,Xb', Q,x m i n ) to depend minimally on them. To 
determine the optimal values fl N and x^ lin , we thus solve the extremality 
conditions dp N {x a ,x b ; fl, x min )/dVl = and dp N (x a ,x b ; ri,x min )/dx min = 0. 
If no extremal point is found, higher derivatives can be used^. The associated 
density matrix pjq is found by normalizing p^: 

/ \ PN \Xa, Xb, ii \X a , Xb, X m [ n yX a , Xb) ), X m - ln yX a , Xb ) ) ... 

PN(X a , X b ) = +00 , TTf — ( N -• I 4 ) 

J_ oo dxp N (x, x; fl N (x, x, x^ in {x, x)),xg in {x, x)) 



x \ )n,ZL = \Pn* 


min V^O.1 •^b)\ 


x b ,%P 

x / Vx]J 
J.. a fc=1 


I dr k V int (x(T k )) 
Jo 
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3 Smearing Formula 

To determine the expectation values ([3]), we proceed as in Ref. LM, and intro- 
duce the Fourier identity 

Tip +00 +00 , -. 

Au*[x]= fdr f dzV int (z) f Je^expU / dr' j(k,r')x(r')\ (5) 



(i 



2n I ^ Jo 



with the current j{k,r') = —ihkS(T — r'). Thus the expectation values ([3]) 
are reduced to path integrals for the euclidean harmonic propagator of ex- 
ponentials of the fluctuating variables. But harmonic expectation values of 
products of exponentials can be written as products of exponentials of pair 
correlation functions as a simple generalization of Wick's rule (the simplest 
well-known application of this being the calculation of the Debye- Waller fac- 
tor for harmonic phonons). Performing the remaining Gaussian integrals we 
obtain directly the smearing formula 



MPnt[*]>n:£„ = II 

fc=l L 

1 A f 1 



h/3 r-\-oo 

dr k dz k V int (z k + x min ) 

J -co 



n exp 1 2 ^ fc ~ Xcl ( Tk ^ G u( T k, T i) [zi - x c i(n)} > ,(6) 



v /(27r)MetG ^ 

where x c \(t) is the classical harmonic path and G denotes the n X n- matrix 
of harmonic Green functions 

^ COShfi(|T fe -Ti| -7i/3) -coshOfo + t ; - ?i/3) 
GHfeTl) -2M sfoh^fi • (?) 

An important advantage of the smearing formula (j3]) over the conventional 
diagrammatic perturbation expansions is that it also allows to calculate cor- 
relation functions for nonpolynomial problems. 



4 Electron Distribution in Hydrogen Atom 

As an example, consider the electron distribution in the hydrogen atom for 
different temperatures shown in Fig.[T] Denoting by r a = \r a \ the distance of 
the electron from the proton, one usually plots the so-called pair distribution 
function gi(r a ), which is equal to (2irf3) 3 ' 2 times the unnormalized density 
Pi(r a ,r a ). Our result is better than an earlier approximation obtained from 
a simple smearing of the effective classical potential of the system ' 
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2.0 4.0 6.0 8.0 10.0 12.0 14.0 



Figure 1. Radial distribution function for an electron in a hydrogen atom. The first-order 
results obtained with isotropic (dashed curves) and anisotropic (solid) variational pertur- 
bation theory are plotted against Storer's resu lts Is l (dotted) and an earlier curve obtained 
from the first-order effective classical potential \±isi (dash-dotted) . 
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